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Abstract 

We discuss the theoretical implications of the scaling properties of 
i*2 at small x observed in recent HERA experiments. We show that 
low Q 2 data display double scaling violations which are adequately 
described by NL InQ 2 corrections. Scaling violations due to summa- 
tions of leading and subleading In - beyond NLO in In Q 2 are however 
disfavoured by the data. 
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The behaviour of the structure function F 2 as a function of x and Q 2 
has been determined [l| by the HI collaboration to an accuracy and with 
a kinematic coverage which makes precision tests of perturbative QCD now 
possible. Whereas the bulk of the data confirm the previously established 
double asymptotic scaling behaviour |2j predicted by perturbative QCD || 
at the leading log level, a study of scaling violations allows one to determine 
the perturbative mechanism which drives the evolution of F 2 as well as to 
disentangle it from its input nonperturbative shape: in particular, it is now 
possible to establish the respective roles of leading logs of Q 2 and - in the 
QCD evolution equations in the HERA kinematic region. 

The contribution of all logarithmic terms of the form 

^(io g QY(iogir. (i) 

to the evolution of structure functions can be included by solution of appro- 
priate renormalization group equations. Double asymptotic scaling follows 
from the symmetric summation of all terms with p = q = r. The recent HI 
data, rescaled by this p rediction,f*| are plotted in fig. 1 versus the product 
of the two scales a = ^/ln ^ In in which \nF 2 should grow linearly with 

a calculable universal slope, and their ratio p = y^n ^ / In ^ , of which F 2 

should be independent (with t = In j^). Whereas the data agree very well 
asymptotically with double scaling, the observed rise of F 2 with a appears to 
be subasymptotically somewhat smaller (i.e. the rescaled F 2 drops); further 
scaling violations are displayed by data with lower values of Q 2 . 

At moderate values of Q 2 , logarithmic contributions to the evolution 
equations with p > q in ( |ID can be important. It turns out that all con- 
tributions with p = r = q + 1 vanish, hence the most important correction 
(corresponding to the inclusion of the most singular contribution in - to the 
NL InQ 2 terms) is given by summing terms with p = q + l = r + l. The 
predicted asymptotic behaviour is still given by a universal double scaling 
form, but receives now a NL (but scheme independent) enhancement, sup- 
pressed by a factor of pa s {Q 2 ) [|J. This correction thus leads to a rise of F 2 
at large p, and a subasymptotic reduction of around 10% in the slope of the 
rise with a in the HERA region. Indeed, if the data are rescaled with this 

*The first subasymptotic correction to double scaling is also included in the rescaling 
factor Rp. 
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Figure 1: The HI data 1 rescaled to the LO double asymptotic scaling prediction. The 
rescaling factor is as in ref. 2. Only data with p > 1, a > 1 and Q 2 > 5 GeV 2 are shown. 
The circles correspond to data with Q 2 = 5 GeV 2 . 



NLO prediction (fig. 2) the predominant double scaling violation apparent 
in fig. 1, namely the drop in a, is removed, and the slope of the rise of 
is now in excellent agreement with QCD (l| . The NLO corrections also have 
the effect of raising somewhat the optimal value of the starting scale. 

The remaining scaling violations displayed in fig. 2 could be simply sub- 
asymptotic corrections, such as contributions to perturbative evolution with 
no In - enhancement (i.e. r = 0), or the influence of the boundary condi- 
tions, or heavy quark thresholds, which can all be kept under control by a 
full two loop analysis. If we wish to use as measured at HERA to perform 
precision tests of QCD it is important to establish whether this is the whole 
story, or whether instead part of the observed double scaling violations may 
be due to higher order corrections in In-, since the treatment of these is 
subject to sizable uncertainties related to scheme dependence |Q- 

The summation of In - contributions can be achieved to all orders [|9] 
by reorganising the perturbative expansion ||, fTj in such a way that In - 
is considered to be leading. It is thus possible to define various expansion 
schemes, each of which will be more accurate in a different kinematic region. 
Whereas in the usual large-x expansion only InQ 2 is leading (i.e. in LO 
p = q > r, in NLO p = q + 1 > r and so on) an expansion more tilted 
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Figure 2: Same as fig. 1, but with NLO rescaling. The rescaling factor is as in ref. 5, 
with Q = 1.8 GeV. The curves have fixed p = 2,3,4 in the a plot and a =1, 1.2, 1.4 in 
the p plot, and correspond to the two loop best fit. 



towards small x (L-expansion) can be constructed by preserving the large 
x form of the LO anomalous dimensions, but adding to the large x NLO 
expressions the higher order leading singularities, i.e. terms with k extra 
powers of a s accompanied by k powers of ln^. It is also possible to define 
an expansion || (double leading or DL expansion) that treats the two logs 
on the same footing: in LO each power of a s is accompanied by either of the 
two logs, while in NLO an overall extra power of a s is allowed.^] 

Such expansion schemes will be adequate provided x is small enough. 
They can then be matched to the large x expansion by imposing continuity 
of anomalous dimensions and coefficient functions at a reference value x = xq. 
This introduces an expansion scheme ambiguity which, once the expansion 



to be used in the small x region has been chosen, is parametrized by £o.[|TO 
If we take for simplicity xq to be Q 2 -independent, this simply means that 
the log which is being summed is actually In a when x < xo, so as xq — > all 
schemes reduce to the standard loop expansion. It should then be possible 
to determine Xq by comparing the computed scaling violations to the data. 



^The NLO anomalous dimensions in this scheme are known explicitly in the quark 
sector [||, but in the gluon sector they can be fixed by choosing factorization schemes in 
which momentum is conserved B. 
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Figure 3: The \ 2 of the three parameter fit (166 d. f.) as a function of x in MS (a) 
and DIS (b) schemes. The solid and dashed lines correspond to the DL expansion with 
standard 6 and Qo 11 factorization, and the dotted line to the L expansion. 



Because of the importance of leading In - effects already at one and two loops 
in the HERA region, one might naively expect the optimal value of Xq to be 
large enough that In - effects beyond two loops are already relevant there. 
Previous data did not allow |10| a determination of Xq, essentially because 
even in the DL scheme the dominant asymptotic behaviour is still given || 
by double scaling. The recent data 0, however, besides being very precise, 
extend well in the subasymptotic (low x and low Q 2 ) region and are thus 
more sensitive to Xq. 

In order to search for higher order logarithms we perform fits to the F 2 
data |l| using two free parameters, namely the exponents X q and X g which 
characterize the small-x behaviour x Xl of the input singlet quark and gluon 
distributions.^] A full set of NLO fits can then be performed with different 
values of xq and in a variety of factorization and small x expansion schemes. 

The dependence of the \ 2 °f these fits on the value of xq is displayed in 
fig. 3. The values of X q , X g and even a s (M z ) are refitted independently at each 
Xq. The quality of the fit always gets monotonically worse as x is increased: 
the data generally favour a very small value of x$. The deterioration of the 
fit is less pronounced in factorization schemes (such as Qo EH compared to 
standard || factorization) or expansion schemes (such as L compared to DL) 
which are closer to the conventional large x two loop expansion. We conclude 



+More details of the fitting procedure are given in ref. 12. 
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Figure 4: The best-fit values of the low- a: exponents X q (crosses) and A g (diamonds) and 
the x 2 01 the ht as a function of the starting scale Qq. 



that the recent data [I] seem to dislike scaling violations which go beyond 
standard two-loop effects. 

We thus focus on the case xq = (two loops). The results of the fit 
as the starting scale is varied, but with a s kept fixed at its best-fit value 
a s (M z ) = 0.122 [12] are shown in fig. 4. The data turn out to require 



soft boundary conditions, with a somewhat softer gluon distribution, and 
a moderately rising quark. As the starting scale Qo is raised, the gluon 
becomes more singular and eventually leads the quark. The quality of the 
fit is largely independent of Qo within a broad range of values. However, 
it deteriorates rapidly if Qo is too high, because the assumed power-like 
behaviour of the input spoils double scaling, or if it is too low, because the 
perturbatively generated rise becomes too strong. The best-fit prediction 
(with Q = 2 GeV) is compared directly to the data in fig. 2. 

We can also look for higher twist contributions by repeating the fits with 
F 2 reparametrized as F^ T = i^ T (l + Cht/Q 2 )- Taking Cut to be x inde- 
pendent its best-fit value turns out to be Cut = 0.2 ± 0.2 GeV 2 . 

In conclusion, our analysis of the 1994 HI data suggests that this is an 
ideal place to perform high-precision tests of QCD: due to the smallness of 
contributions related to higher logs of - and higher twists, the usual NLO 
perturbative expansion is perfectly adequate. The absence of scaling viola- 
tions related to higher logs of as expressed in the unnaturally small value 
of xo, is however unexpected, suggesting that our understanding of the way 
these logs should be summed is as yet incomplete. 
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